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ABSTRACT 

We investigate observational constraints on the cosmic equation of state from 
measurements of angular size for a large sample of milliarcsecond compact radio- 
sources. The results are based on a flat Friedmann- Robert son- Walker (FRW) 
type models driven by non-relativistic matter plus a smooth dark energy compo- 
nent parametrized by its equation of state p x = up x (—1 < uj < 0). The allowed 
intervals for u and Q m are heavily dependent on the value of the mean projected 
linear size I. For I ~ 20h~ 1 — 30/i _1 pc, we find fi m < 0.62, uj < —0.2 and 
Q m < 0.17, u < —0.65 (68% c.l.), respectively. As a general result, this analysis 
shows that if one minimizes \ 2 f° r the parameters I, Q m and uj, the conventional 
flat ACDM model (u = -1) with fi m = 0.2 and I = 22.6/r 1 pc is the best fit for 
these angular size data. 

Subject headings: cosmology: theory - dark matter - distance scale 

A large number of recent observational evidences strongly suggest that we live in a flat, 
accelerating Universe composed by ~ 1/3 of matter (barionic + dark) and ~ 2/3 of an exotic 
component with large negative pressure, usually named dark energy or "quintessence" . The 
basic set of experiments includes: observations from SNe la (Perlmutter et al. 1998; 1999; 
Riess et al. 1998), CMB anisotropics (de Bernardis et al. 2000), large scale structure (Bah- 
call 2000), age estimates of globular clusters (Carretta et al. 2000; Krauss 2000; Rengel et al. 
20001) and old high redshift galaxies (OHRG's) (Dunlop 1996; Krauss 1997; Alcaniz & Lima 
1999; Alcaniz & Lima 2001). It is now believed that such results provide the remaining piece 
of information connecting the inflationary flatness prediction (fix = 1) with astronomical 
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observations, and, perhaps more important from a theoretical viewpoint, they have stimu- 
lated the current interest for more general models containing an extra component describing 
this unknown dark energy, and simultaneously accounting for the present accelerated stage 
of the Universe. 

The absence of a convincing evidence concerning the nature of this dark component 
gave origin to an intense debate and mainly to many theoretical speculations in the last few 
years. Some possible candidates for "quintessence" are: a vacuum decaying energy density, 
or a time varying A-term (Ozer & Taha 1987; Freese 1987; Carvalho et al 1992, Lima and 
Maia 1994), a relic scalar field (Peebles & Ratra 1988; Frieman et al 1995; Caldwell et al 
1998; Saini et al 2000) or still an extra component, the so-called "X-matter", which is simply 
characterized by an equation of state p x = ujp x , where uj > — 1 (Turner & White 1997; Chiba 
et al 1997) and includes, as a particular case, models with a cosmological constant (ACDM) 
(Peebles 1984). For "X-matter" models, several results suggest Q x ~ 0.7 and uj < —0.6. 
For example, studies from gravitational lensing + SNe la provide uj < —0.7 at 68% c.l. 
(Waga & Miceli 1999; see also Dev et al. 2001). Limits from age estimates of old galaxies 
at high redshifts require uj < —0.27 for Q m ~ 0.3 (Lima & Alcaniz 2000a). In addition, 
constraints from large scale structure (LSS) and cosmic microwave background anisotropies 
(CMB) complemented by the SN la data, require 0.6 < Q x < 0.7 and uj < —0.6 (95% c.l.) 
for a flat universe (Garnavich et al 1998; Perlmutter et al 1999; Efstathiou 1999), while for 
universes with arbitrary spatial curvature these data provide uj < —OA (Efstathiou 1999). 

On the other hand, although carefully investigated in many of their theoretical and ob- 
servational aspects, an overview on the literature shows that a quantitative analysis on the 
influence of a "quintessence" component {uj = p x /p x ) in some kinematic tests like angular 
size-redshift relation still remains to be analysed. Recently, Lima & Alcaniz (2000b) studied 
some qualitative aspects of this test in the context of such models, with particular emphasis 
for the critical redshift z m at which the angular size takes its minimal value. As a general 
conclusion, it was shown that this critical redshift cannot discriminate between world mod- 
els since different scenarios can provide similar values of z m (see also Krauss & Schramm 
1993). This situation is not improved even if evolutionary effects are taken into account. In 
particular, for the observationally favoured open universe (f2 m = 0.3) we found z m = 1.89, 
a value that can also be obtained for quintessence models having 0.85 < £l x < 0.93 and 
— 1 < uj x < —0.5. Qualitatively, it was also argued that if the predicted z m is combinated 
with other tests, some interesting cosmological constraints can be obtained. 

In this letter, we focus our attention on a more quantitative analysis. We consider the 
9 — z data of compact radio sources recently updated and extended by Gurvits et al. (1999) 
to constrain the cosmic equation of state. We show that a good agreement between theory 
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and observation is possible if Q m < 0.62, uj < —0.2 and £l m < 0.17, uj < —0.65 (68% c.l.) 
for values of the mean projected linear size between I ~ 20/i -1 — 30k" 1 pc, respectively. In 
particular we find that a conventional cosmological constant model (u = —1) with Q m = 0.2 
and I = 22.64/i~ 1 pc is the best fit model for these data with x 2 — 4.51 for 9 degrees of 
freedom. 

For spatially flat, homogeneous, and isotropic cosmologies driven by nonrelativistic mat- 
ter plus an exotic component with equation of state, p x = ujp x , the Einstein field equations 
take the following form: 



I )2 = H ° 



fi m (§) 3 + ^(^) 3(1+w) 



R = ~2 H ° 



fi m (§) 3 + (3c + iR(§) 3(1+w) 



:i) 



(2) 



where an overdot denotes derivative with respect to time, H Q = 100/iKms~ 1 Mpc~ 1 is the 
present value of the Hubble parameter, and fl m and fl x are the present day matter and 
quintessence density parameters. As one may check from (1) and (2), the case cu — — 1 
corresponds effectively to a cosmological constant. 

In such a background, the angular size-redshift relation for a rod of intrinsic length / 
can be written as (Sandage 1988) 

In the above expression D is the angular-size scale expressed in milliarcseconds (marcs) 

D = ™* (4) 

c 

where / is measured in parsecs (for compact radio-sources), and the dimensionless coordinate 
£ is given by (Lima & Alcaniz 2000b) 

/ r- (5) 

J(i+z)-^ x [n^- 1 + (i - n ra )x-( 1 + 3 ^] * 



The above equations imply that for given values of I, fi m and u, the predicted value of 
6(z) is completely determined. Two points, however, should be stressed before discussing 
the resulting diagrams. First of all, the determination of fl m and uj are strongly dependent 
on the adopted value of I. In this case, instead of assuming a especific value for the mean 
projected linear size, we have worked on the interval / ~ 20/i _1 — 30/iT 1 pc, i.e., I ~ O(40) 
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pc for h = 0.65, or equivalently, D = 1.4 — 2.0 marcs. Second, following Kellermann (1993), 
we assume that possible evolutionary effects can be removed out from this sample because 
compact radio jets are (i) typically less than a hundred parsecs in extent, and, therefore, 
their morphology and kinematics do not depend considerably on the intergalactic medium 
and (ii) they have typical ages of some tens of years, i.e., they are very young compared to 
the age of the Universe. 

In our analysis we consider the angular size data for milliarcsecond radio-sources recently 
compiled by Gurvits et al. (1999). This data set, originally composed by 330 sources 
distributed over a wide range of redshifts (0.011 < z < 4.72), was reduced to 145 sources 
with spectral index —0.38 < a < 0.18 and total luminosity Lh 2 > 10 26 W/Hz in order 
to minimize any possible dependence of angular size on spectral index and/or linear size on 
luminosity. This new sample was distributed into 12 bins with 12-13 sources per bin (see Fig. 
1). In order to determine the cosmological parameters Q m and uj, we use a x 2 minimization 
neglecting the unphysical region f2 m < 0, 



where 6(zi, /, Q m , uj) is given by Eqs. (3) and (5) and 9 oi is the observed values of the angular 
size with errors <7j of the ith bin in the sample. 

Figure 1 displays the binned data of the median angular size plotted against redshift. 
The curves represent flat quintessence models with Q m = 0.3 and some selected values of 
uj. As discussed in Lima & Alcaniz (2000b), the standard open model (thick line) may be 
interpreted as an intermediary case between ACDM (uj = —1) and quintessence models with 
uj < —0.5. In Fig. 2 we show contours of constant likelihood (95% and 68%) in the plane 
uj-Q m for the interval I ~ 20k- 1 -SOh' 1 pc. For I = 20.58/i" 1 pc (D = 1.4 marcs), the best 
fit occurs for Q m = 0.26 and uj = —0.86. As can be seen there, this assumption provides 
Q m < 0.48 and uj = —0.3 at la. In the subsequent panels of the same figure similar analyses 
are displayed for / ~ 22.05/i _1 pc (D = 1.5 marcs), I ~ 23.53/i -1 pc (D = 1.6 marcs) and 
I ~ 29.41/i -1 pc (D = 2.0 marcs), respectively. As should be physically expected, the limits 
are now much more restrictive than in the previous case because for the same values of 9 oi 
it is needed larger £(z) (for larger /) and, therefore, smaller values of uj. For / ~ 29.41/iT 1 
pc, we find Q m = 0.04 and uj = — 1 as the best fit model. For intermediate values of 
namely, / = 22.0/i" 1 pc (D = 1.5 marcs) and / = 23.5h~ 1 pc (D = 1.6 marcs), we have 
f2 m = 0.22, uj = —0.98 and f2 m = 0.16 and uj = —1, respectively. In particular, for smaller 
values of /, e.g., / ~ 14.r0/i~ 1 pc (D = 1.0 marcs) we find Q m = 0.36, uj = -0.04. As 
a general result (independent of the choice of /), if we minimize x 2 f° r U an d uj, we 
obtain / = 22.64/r 1 pc (D = 1.54 marcs), fi m = 0.2 and u = -1 with x 2 = 4.51 for 9 
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degrees of freedom (see Table 1). It is worth notice that our results are rather different 
from those presented by Jackson & Dodgson (1996) based on the original Kellermann's data 
(Kellermann 1993). They argued that the Kellermann's compilation favours open and highly 
decelerating models with negative cosmological constant. Later on, they considered a bigger 
sample of 256 sources selected from the compilation of Gurvits (1994) and concluded that the 
standard flat CDM model is ruled out at 98.5% confidence level whereas low-density models 
with a cosmological constant of either sign are favoured (Jackson & Dodgson 1997). More 
recently, Vishwakarma (2001) used the updated data of Gurvits et al. (1999) to compare 
varying and constant ACDM models. He concluded that flat ACDM models with £7 m = 0.2 
are favoured. 

At this point it is also interesting to compare our results with some recent determina- 
tions of uj derived from independent methods. Recently, Garnavich et al. (1998) using the 
SNe la data from the High-Z Supernova Search Team found uj < —0.55 (95% c.l.) for flat 
models whatever the value of Q m whereas for arbitrary geometry they obtained uj < —0.6 
(95% c.l.). As commented there, these values are inconsistent with an unknown component 
like topological defects (domain walls, string, and textures) for which uj — — |, being n the 
dimension of the defect. The results by Garnavich et al. (1998) agree with the constraints 
obtained from a wide variety of different phenomena (Wang et al. 1999), using the "concor- 
dance cosmic" method. Their combined maximum likelihood analysis suggests uj < —0.6, 
which is less stringent than the upper limits derived here for values of I > 20h~ 1 pc. More 
recently, Balbi et al. (2001) investigated CMB anisotropies in quintessence models by us- 
ing the MAXIMA-1 and BOOMERANG-98 published bandpowers in combination with the 
COBE/DMR results (see also Corasaniti & Copeland 2001). Their analysis sugests Q x > 0.7 
and — 1 < uj < —0.5 whereas Jain et al (2001) found, by using image separation distribution 
function of lensed quasars, —0.75 < uj < —0.42, for the observed range of Q m ~ 0.2 — 0.4 
(Dekel et al. 1997). These and other recents results are summarized in Table 2. 

Let us now discuss briefly these angular size constraints whether the adopted X-matter 
model is replaced by a scalar field motivated cosmology, as for instance, that one proposed by 
Peebles and Ratra (1988). These models are defined by power law potentials, V((f>) ~ <p~ a , 
in such a way that the parameter of the effective equation of state = p<f>/p<f>) may become 
constant at late times (or for a given era). In this shown elsewhere (Lima & Alcaniz 

2000c), the dimensionless quantity £ defining the angular size reads 



Comparing the above expression with (5) we see that if u — —2/(2 + a) this class of models 
may reproduce faithfully the X-matter constraints based on the angular size observations 
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presented here. However, as happens with the Supernovae type la data (Podariu & Ratra 
2000), the two sets of confidence contours may differ significantly if one goes beyond the 
time independent equation of state approximation. Naturally, a similar behavior is expected 
if generic scalar field potentials are considered. 

Finally, we stress that measurements of angular size from distant sources provide an 
important test for world models. However, in order to improve the results a statistical 
study describing the intrinsic lenght distribution of the sources seems to be of fundamen- 
tal importance. On the other hand, in the absence of such analysis but living in the era 
of precision cosmology, one may argue that reasonable values for astrophysical quantities 
(like the characteristic linear size /) can be infered from the best cosmological model. As 
observed by Gurvits (1994), such an estimative could be useful for any kind of study en- 
volving physical parameters of active galactic nuclei (AGN). In principle, by knowing / and 
assuming a physical model for AGN, a new method to estimate the Hubble parameter could 
be established. 
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Fig. 1. — Angular size versus redshift for 145 sources binned into 12 bins (Gurvits et al. 
1999). The curves correspond to the characteristic linear size / = 22.64/i _1 pc. Thick solid 
curve is the prediction of the standard open model (fi m = 0.3). 
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Matter density parameter (^ m ) 

Fig. 2. — Confidence regions in the u — Q m plane according to the updated sample of angular 
size data (Gurvits et al. 1999). The solid lines in each panel show the 95% and 68% likelihood 
contours for flat quintessence models. 
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D (mas) Ih (pc) f2 m uo \ 2 

L4 20.58 0.26 -0.86 4.56 

1.5 22.05 0.22 -0.98 4.52 

1.6 23.53 0.16 -1 4.54 
2.0 29.41 0.04 -1 5.57 

Best fit: 1.54 22.64 0.2 -1 4.51 



Table 1: Limits on uj from 9 — z relation 



Method Author Q m u 

CMB+SNe la.. Turner & White (1997) ~ 0.3 ~ -0.6 

Efstathiou (1999) ~ < -0.6 

SNe la Garnavich et al. (1998) ~ < -0.55 

SGL+SNe la.. Waga & Miceli (1999) 0.24 < -0.7 

SNe Ia+LSS... Perlmutter et al. (1999) ~ < -0.6 

Various Wang et al. (1999) 0.2 - 0.5 < -0.6 

OHRG's Lima & Alcaniz (2000a) 0.3 < -0.27 

CMB Balbi et al. (2001) 0.3 < -0.5 

Corasaniti & Copeland (2001) ~ < -0.96 

SGL Jain et al. (2001) 0.2 - 0.4 > -0.75, < -0.55 



Table 2: Limits to u for a given Q 
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